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Abstract 



O 

We establish a maximin characterisation of the linear escape rate of the orbits of 
a non-expansive mapping on a complete (hemi-)metric space, under a mild form of 
i— i Busemann's non-positive curvature condition (we require a distinguished family of 

geodesies with a common origin to satisfy a convexity inequality) . This characterisa- 
tion, which involves horof unctions, generalises the Collatz-Wielandt characterisation 
of the spectral radius of a non-negative matrix. It yields as corollaries a theorem of 
Kohlberg and Neyman (1981), concerning non-expansive maps in Banach spaces, a 
variant of a Denjoy- Wolff type theorem of Karlsson (2001), together with a refine- 
ment of a theorem of Gunawardena and Walsh (2003), concerning order-preserving 
positively homogeneous self-maps of symmetric cones. An application to zero-sum 
stochastic games is also given. 

in 

^ 1 Introduction 

A self-map T of a metric space (X, d) is non-expansive if d(T(x), T(y)) ^ d(x, y). A general 
problem consists in studying the asymptotic behaviour of the orbits of T. This is moti- 
vated in particular by the celebrated theorem of Denjoy and Wolff Den26[fvVol26a Wol26b| 



on the iteration of holomorphic self-maps of the unit disk (these maps are non-expansive in 
^ Poincare's metric). Other motivations arise from the cases of non-expansive self-maps of 

Banach spaces and of self-maps of cones that are non-expansive in Hilbert's, Thompson's 
or Riemannian metric, which have received a considerable attention, due in particular to 



applications in game theory |RS01,NS04 , discrete event systems Gun03|, quadratic opti- 



mal control or filtering Bou95|, and non-linear Perron- Frobenius theory Nus88 . Several 
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Denjoy- Wolff type results, either in the setting of metric spaces, or concerning the special 
case of cones have appeared, see in particular Bea97,Kar01 , AGLN06 Nus07 Lin07 . The 



reader may consult the monograph RS05 for an overview of the field. 



In this paper, we consider the linear escape rate 



p{T) = lim 

k— »oo 



d{x,T k {x)) 



k 



The latter always exists, by a standard subadditive argument. Our main result is the 
following maximin characterisation of the escape rate. 

Theorem 1. Let T be a non-expansive self-map of a complete metrically star shaped 
hemi-metric space (X, d) . Then, 



inf d(y,T(y)) = p(T) 



max inf h(T(x)) — h(x) 

h x£X 



where the maximum, which is attained, is taken over the set of Martin functions h of 
(X,d). If in addition the hemi-metric d is bounded from below (in particular, if d is a 
metric), and if p(T) > 0, then any function attaining the maximum is a horofunction. 

Let us explain the terminology. Hemi-metrics are analogous to metrics, but d(x, y) 
and d(y,x) may differ. Their definition is a variation of that of weak metrics in |PT09| . 
The notion of metrically star shaped space is the object of Definition [5] below. It requires 
a distinguished family of geodesies with a common origin to satisfy a convexity inequality. 



Hence, it is a mild form of Busemann's classical non-positive curvature condition Pap05 



Martin functions and horofunctions are special Lipschitz functions of constant 1 which 
arise when compactifying the hemi-metric space (X, d). The set of horofunctions (the 
horoboundary) provides an abstract boundary of X. See £ |2.3 for details. 

Theorem [T] is inspired by the classical Collatz-Wielandt characterisation of the Perron 
root p(M) of a n x n non-negative matrix M. The latter shows that 



inf max 

yeintR™ lsgi^n 



(My), . (Mu) 
= p\M) = max mm 

yi ueR r l IsSisSn m 



1.2) 



Here, M + denotes the set of real non-negative numbers, and so W]_ is the standard positive 



cone. Nussbaum showed in Nus86 that the latter formula holds more generally when the 



map y h- >■ My is replaced by a non-linear continuous order-preserving T of the standard 
positive cone. When T preserves the interior of this cone, it is non-expansive in the 



reverse) Funk hemi-metric PT09 , Wal08] defined on this interior by 



RFunk(x,y) := log sup — 



Thus, 



RFunk(T(x),T(y)) ^ RFunk(x,?/) 



Vx, y G int; 



The Collatz-Wielandt formula and its non-linear extension in Nus86 will be shown to be 



special instances of the maximin formula (1.1) 
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Theorem[T]is also motivated by the case of a Banach space X. Then, one may consider 
the limit 



lim T k {x)/k 



1.3) 



For instance, in the setting of zero-sum games, T is the dynamic programming or Shapley 



operator Sha53| (which acts on a Banach space of continuous functions equipped with 
the sup-norm), and the latter expression represents the limit of the mean payoff per time 
unit as a function of the initial position, when the horizon of the game tends to infinity, 



see £3.3 for more background. 



A theorem of Pazy Paz71| shows that the limit (1.3) does exist if X is a Hilbert 
space; further improvements, under suitable strict convexity assumptions were done by 
Reich|Rei73] and then by Kohlberg and Neyman KN81]: the limit exists in the weak 



(resp. strong) topology if X is reflexive and strictly convex (resp. if the norm of the dual 
space X* is Frechet different iable). Without strict convexity, the limit may not exist, even 



when X is finite dimensional, however, a general result of Kohlberg and Neyman KN81 



(Corollary 20 below) shows that there is always a linear form (f £ X* of norm 1 such that 



V {T\x)) > p(x) + k-p{T) 



holds for all k £ N, where 



P(T) :- 



inf d(y,T(y)) 

yex 



is precisely the term at the left-hand side of (1.1). The theorem of Kohlberg and Neyman 



may also be thought of as a special case of Theorem [TJ Indeed, any horofunction h 
attaining the maximum in (1.1) satisfies 



h(T k (x)) > h(x) + kp(T) ,WxeX 



for all k £ N. We shall see that the Kohlberg-Neyman theorem follows readily from 
this result. We shall also recover as corollaries a generalisation of a result of Gaubert and 
Gunawardena |GG04] (Corollary [23)) concerning Shapley operators, as well as a refinement 
of a result of Gunawardena and Walsh GW03] valid in symmetric cones, Corollary 32 
below. 



Theorem[T]turns out to be related to a Denjoy- Wolff type theorem of Karlsson KarOl 
Theorem [25] below. Karlsson showed that an inequality similar to (1.4) holds. He used 
a general subadditivity argument, which does not require any non-positive curvature 
condition. The statement of |Kar01| assumes the metric space to be proper, i.e., closed 
balls to be compact, but this assumption can be relaxed by defining horofunctions with 
respect to the topology of pointwise convergence as we do here. However, there are 
essential discrepancies between (1.4) and the result of KarOl : the term p(T) is replaced 



there by the linear escape rate p(T), and the horofunction h depends on the choice of the 
point x, whereas it does not in ( 1.4[ ). We shall see that without a non-positive curvature 
condition, the formula (1.4) may not hold (we give an example in which p(T) > p(T) 
and h necessarily depends on x). Hence, the maximin Theorem [I] holds only under more 
restrictive circumstances than the Denjoy- Wolff type result of Karlsson. By comparison, 
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the interest of the more special Theorem [T] lies in its strong duality nature: it allows one 
to certify that p(T) < a, or dually, that p(T) > f3, by exhibiting an element y G X such 
that d(y,T(y)) < a, or dually a Martin function h such that mf xe x h(T(x)) — h(x) > (3. 
It should be noted in this respect that when X is a Banach space with a polyhedral norm, 
the horoboundary admits a simple effective description, see KMN06 , Wal07 , AGW09| . 



2 Definitions and preliminary results 
2.1 Metrically star-shaped spaces 

Definition 2. We say that 6 : X xX — )■ R is a hemi-metric on a set X if the two following 
conditions are satisfied for all (x,y, z) G X 3 : 

a) 8{x,z) ^ 8(x,y) + 5(y,z) 

b) 5(x, y) = 5(y, x) — if and only if x = y. 
We then say that (X, 5) is a hemi-metric space. 

Notice that a hemi-metric is generally not a metric, since it is neither symmetric nor 



m 



non-negative. This definition is closely related to the one of a weak metric in PT09 
which 5 is required to be non-negative. The last condition of the definition corresponds to 
the weak separation condition in the latter reference. We allow 5 to take negative values in 
order to deal with order-preserving positively homogeneous self-maps of cones (Section 
To any hemi-metric, one can canonically associate a metric by the following lemma. 

Lemma 3. For any hemi-metric 5, the function d(x, y) = max(S(x, y), 5(y, x)) is a metric 
on X. 

Proof. Verifying that d is symmetric and satisfies the triangular inequality is easy. The 
positivity of d comes from the fact that for all x and y, = 5(x,x) ^ S(x,y) + 5(y,x). 
Finally, if d(x,y) = then both S(x,y) and S(y,x) are non-positive so by the same 
argument they are both null, thus x = y. □ 

In the sequel, X is equipped with the topology induced by the metric d. We shall say 
that (X, 5) is complete when the associated metric space (X, d) is complete. 

Definition 4. A geodesic joining a point x G X to a point y G X is a map 7 : [0, 1] — > X 
such that 7(0) = x, 7(1) = y, and such that for all ^ s ^ t ^ 1, 

6(7(s),y(t)) = (t-s)6(x,y). 

Definition 5. We say that (X, 5) is metrically star-shaped with centre x° if there exists 
a family of geodesies {7 J/ }j /g x, such that y y joins the centre x° to the point y, and such 
that the following inequality is satisfied for every (y,z) G X 2 and s G [0,1]: 

5( ly (s), lz (s)) ^ s5(y,z). (2.1) 

If the hemi-metric 5 is not a metric, we also require that for any y, the quantity 8(y, 7j,(s)) 
tends to as s goes to 1. 
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The condition (2.1 ) is a form of metric convexity Pap05| . In particular any Busemann 
space Pap05| is metrically star-shaped, but our definition is less demanding since we only 
require the inequality (2.1) to be satisfied for one specific choice of geodesies. 



Example 6. Any Banach space (X, || ■ ||) is metrically star-shaped with respect to any 
centre x° : it suffices to take the straight lines as geodesies, i.e., for any choice of centre 
x°, the choice of j y (s) = x° + s(y — x°) yields a metrically star shaped space. Notice that 



some Banach spaces are not Busemann spaces Pap05 



The next example relies on the following notion. 

Definition 7. A map p from a Banach space (X, \\ 
written as 



p{z) 



sup <p{z) 



to IR is a hemi-norm if it can be 



(2.2) 



where E is a bounded subset of the dual space X*, and 

p(z) = p(—z) = =>- 







(2.3) 



A hemi-norm is always Lipschitz. It is easily verified that if p is a hemi-norm, then 
S(x,y) := p(y — x) is a hemi-metric. Note that Condition (2.3) holds if and only if the 
orthogonal set of E is reduced to the zero vector. As in the case of Banach spaces, 
(X, 5) is metrically star shaped: it suffices to choose the straight lines as geodesies. We 
shall say that the hemi-norm is compatible with the Banach space (X, || ■ ||) if ||,z|| = 
ma.x(p(z) , p(— z)). Then, the metric d obtained from 5 is the one associated to the norm 
of X, and in particular, the hemi-metric space (X, 5) is complete. 

Example 8. The norm of the Banach space X is a special case of compatible hemi-norm, 
in which 

p(z) = \\z\\ = sup ip(z) 



where B* is the dual unit ball of X*. Then, by Bauer maximum principle AB06, The- 



orem 7.69], the supremum is attained at some extreme point of B* , so the set E arising 



in the definition of the hemi-norm (see Equation 2.2) can be taken to be either B* or the 
set of extreme points of B*. 

Example 9. A useful example of hemi-norm is the following. Let Q denote a compact 
topological space, let X = ^(Q) denote the Banach space of continuous functions from Q 
to R, equipped with the sup-norm. Consider 5(x, y) = t(y — x), where t denotes the "top" 
operator which gives the maximum of a function, i.e., 



t(x) 



maxa;(w) 



The metric d obtained from 5 is the sup-norm. We can take for E the set of evaluations 
functions {e w ,w G fl} with e LU (x) = x(oj). 

Example 10. Another class of examples concern symmetric cones, which include cones of 
positive semidefinite matrices. We shall discuss them further in Section |4} 
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2.2 Non-expansive mappings 



In this section we consider a metrically star-shaped space (X, S) with centre x°, as well 
as a map T : X — > X that is non-expansive with respect to the hemi-metric 5, meaning 
that for all (x,y) G X 2 , 

5(T(x),r(y))0(z,y). 

Definition 11. To each non-expansive mapping T, we associate the two following quan- 
tities: 

p{T) = m£6(x,T{x)) , (2.4) 

, s 5(x,T k (x)) 5(x,T k (x)) , s 

o T = lim V - , = inf V ' , ;; . 2.5 

k^+oo k k^l fc 

Thus, the number p(T) measures the linear escape rate of the orbits of T. It is well 
defined. Indeed, for any x G X, since T is non-expansive, the sequence = 
is su badditiv e, meaning that Uk+i ^ + Ui for all k, I ^ 1, and so, a classical argument 
Lemma 1.18]) shows that the limit lim^ +OC) fc -1 <5 (a;,T fc (x)) does exist and 



Bow08 



[see 

is equal to inf&^i k~ l 5 (x, T k (x)). Moreover, since 

5{x,T k (x)) < tffoy) +tf(y,T fc (i/))+ 5 (T fc (y),T fc (x)) 
< 5(x,y) + 5(y,x) + 5(y,T k (y)) 

the limit is independent of the choice of x G X. 

Lemma 12. 27ie following inequality is satisfied for any non-expansive mapping T: 

P(T) < ~p{T). 

Proof. Observe that for any x G X and k ^ 1, 

fe-i 

5 (x,T fc (a;)) ^ 5 (T^a;), T /+1 (x)) ^ M(x,T(x)) 
and take the infimum on both k and x. □ 



2.3 The horofunction boundary 



The horofunction boundary of a metric space was defined by Gromov |Gro81 . See 



also BGS85 , Bal95, Ch. II] and Rie02 . In fact, the same construction can be performed 
with a hemi-metric. The details can be found in | A GW09] , in which a horofunction-like 
boundary is defined for discrete optimal control problems, by analogy with the Martin 
compactification arising in the theory of Markov processes. See also |IM07|. 

Let us fix an arbitrary point x G X (the basepoint). We define a map i from X to the 
set functions from X to M by associating to any x G X the following function i(x): 



i(x) : y [i(x)](y) := S(x,x) - 6(y,x). 



6 



Using the triangular inequality, we get that, for all x,y G X, 



-S(y,x) ^ [i{x)]{y) < 6{x,y) 

Hence, the set i(X) := {i(x) \ x G X} can be identified to a subset of the product space 
rij/ext - fi{y> <K^> 2/)] Tychonoff's theorem, the later space is compact for the 

product topology. Hence, the closure of i(X) in the topology of pointwise convergence 
(which is the same as the product topology) is compact. We denote by jtft this closure. 
The elements of the boundary J%? := \ i(X) are called horof unctions. We will call 
Martin functions the elements of jtft (the Martin space). Note that the choice of the 
basepoint is irrelevant (changing the basepoint translates all the Martin functions by the 
same constant). 

Remark 13. The term of horofunction is sometimes used in a more restrictive sense, to 
denote what is also called a Busemann function, which is a horofunction obtained as 
the limit of a family of functions (i(x s )) s ^o taken along an "infinite geodesic" (x s ) s ^q, 



see |Rie02 . Note also the sign in the definition of i(x) (which implies, in the Banach 
space case, that the map i(x) is concave). This sign is chosen consistently with potential 
theory (the opposite of the distance is the analogue of the Martin kernel). In the metric 
geometry literature, the opposite choice is most often made. 

Remark 14. The horofunctions are often defined as the closure of i(X) in the topology 
of uniform convergence on bounded sets; indeed, the injection x i— > i(x) is continuous for 



this topology, and this injection is an embedding if X is a complete geodesic space Bal95 
Ch. II]. It can be checked that every map i(x) is Lipschitz of constant 1 with respect 
to the metric d, and so, by the Ascoli-Arzela theorem, the closure of i(X) in this sense 
is compact as soon as X is proper (meaning that every closed ball is compact). In the 
present work, we do not require X to be proper, but define rather as the closure of i(X) 
in the topology of pointwise convergence, so that jtft is always compact. Note however 
that the injection x h-> i(x) may not be an embedding from X to ^# (the topology on ^# 
is too weak in general, so the inverse of the map x H- i(x) may not be continuous). 

Remark 15. The topology of the Martin space is metrisable as soon as X is a countable 
union of compact sets, see for instance |AGW09, Remark 7.10]. Then, every horofunction 



is the limit of a sequence of functions i(x m ), where x m is a sequence of elements of X. 



3 The main result and some of its consequences 
3.1 The main result 

We shall derive Theorem [T] from the following result. 

Theorem 16. Let (X, 5) be a complete metrically star-shaped hemi-metric space, and let 
T : X — > X be non- expansive. Then there exists a Martin function h G ^ such that for 
all x G X , 

h(T(x)) ^ h(x)+p(T) . 

Moreover, if 5 is bounded from below (in particular, if 5 is a metric) and ifp(T) > 0, then 
h is necessarily a horofunction. 
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Proof. Let x° and {'y y }y£X as in Definition |5j for any a G [0, 1[ denote by r a : X — > X the 
function y — > 7 y (a) and recall that r a is a-contracting by definition. By completeness, 
we can thus define for any a the point y a G X as the only solution of the fixed point 
equation 

T{r a {y a )) = y a . 

Then for any 

5(x,y a ) - 5 (T(x),y a ) = 5(x,y a ) - 5 (T(x),T(r a (y a ))) 

^ 5(x, y a ) — 5 (x, r a {y a )) (by non-expansiveness) 
^ 5(x, y a ) - 5 (x, r a (x)) - 5 (r a (x),r a (y a )) 
^ (1 — a)5(x, y a ) — 5 (x, r a (x)) (since r a is an a-contraction) 
^ (1 — a)5(x°, y a ) — (1 — a)5(x°, x) — 5 (x, r a (x)) 
= S(r a (y a ),y a ) - (1 - a)S(x°,x) - 5 (x,r a (x)) 
> 5 (T(r a (y a ),T(y a ))) - (1 - a)5(x°, x) - 5 (x, r a (x)) 

(by non-expansiveness) 
= 5 (y a ,T(y a ))) - (1 - a)S(x°,x) - S (x,r a (x)) 
^ p(T)-(l-a)5(x°,x)-5(x,r a (x)). 

Since the space ^ is compact, the family of functions {i{y a ))o<a<i admits a limit point 
h G ^ as a tends to 1. Passing to the limit in the previous inequality, and using the 
additional assumption in Definition |5j we deduce that —h(x) + h(T(x)) ^ p(T). 

Assume now by contradiction that 5 is bounded from below, that p(T) > 0, and 
that h = i(z) for some z G X. Then, we deduce from h(T(x)) ^ h(x) + p(T) and 
h(y) = —S(y, z) + 5(b, z) that 

-5(T k (x),z) + 6(x,z) ^ kp(T) . 

The right-hand side of this expression tends to +oo as k tends to infinity, but the left-hand 
side of this expression is bounded above independently of k, since 5(T h (x), z) is bounded 
from below. This is impossible, hence, h G Jti? is a horofunction. (More generally, the 
same argument shows that for any net of functions (i(x ( a)) / 3 S B converging to h, the net 
(x/3)/3 £ b cannot have a bounded subnet.) □ 

We obtain as immediate consequences the following corollaries: 



Corollary 17. Let X,5 andT be as in Theorem 16 Then, there exists a Martin function 
h G such that for all x G X , 

h(T k (x)) ^ h{x) + kp{T) . (3.1) 

Moreover, if 5 is bounded from below and ifp(T) > 0, then h is necessarily a horofunction. 



Corollary 18. // the assumptions of Theorem 16 are satisfied, then p(T) = p(T). 
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Proof. We deduce from (3.1) that 



m < \{h{T\x))-h{ X )) ^ T * {x)) 



which yields p(T) ^ p(T). Then, the result follows by Lemma 12 



□ 



Theorem [T] follows readily by combining Theorem 16 and Corollary 18 



3.2 The Kohlberg-Neyman theorem revisited 



We now apply Theorem [16] to non-expansive self-maps of a Banach space X, and more 
generally, to those maps that are non-expansive in a compatible hemi-norm, as defined in 
Section 12.11 

Theorem 19. Let (X, || • ||) be a Banach space, and let 

p(z) = sup <fi(z) 

c/>£E 

be a compatible hemi-norm, where E is a bounded subset of the dual space X* . Let T be 
non-expansive for 5. Then, for every iGl, there exists a linear form (fi in the weak-star 
closure of E such that for every k G N, 

<fi(T k (x)) ^<fi(x) + kp(T). 

Proof. Let cloE denote the closure of E in the weak-star topology. Observe that, for all 

u e x, 

p{u) = sup <p(u) . 
</>edo E 

Since E is bounded, cloi? is weak-star compact, and so, the map (fi h-> 4>{u) attains its 
maximum on clo E. 

Fix now x and y m. X. It follows from the previous discussion that for all z G X, 
%)W - Kv)( x ) = P(V -x)- p(y - z) = <p(y -x)- p(y - z) 
for some G clo E, independent of z, and so 

%)W - Kv)( x ) ^ <P(v ~ x )- <i>{y - z ) = <t>{ z - x )- 

In other words, is in the super-differential of i(y) at point x. 



Let h be the Martin function that appears in Corollary 17 , h = limn i(yd) for some net 
{Hd)d&D- By the previous observation, for every d G D we can find an element (fid G cloi? 
such that for every k G N, 

i{y d ){T k {x)) - i{y d ){x) ^ d (T fc (x) - x). 

Taking the limit along some subnet, we can find a linear form (fi in c\oE such that for 
every fceN, 

h(T k (x)) - h(x) «C (fi(T k (x) - x). 



Since by Corollary 17 h(T k (x)) — h(x) ^ kp(T), the theorem is established. □ 
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We now recover the following result. 



Corollary 20 (Kohlberg and Neyman KN81, Ney03| ). Let (X, || ■ ||) be a Banach space, 
let T : X —J- X be non- expansive and assume that p(T) > 0. Then for any x G X , there 
exists a continuous linear form of norm one, such that 



4>(T k {x)) ^ <f>(x) + kp(T) 



(3.2) 



for all fceN. Moreover, can be taken in the weak-star closure of the set of extreme 
points of the dual unit ball. 



Proof. We apply Theorem 19, taking E to be the set of extreme points of the dual unit 
ball. This gives a in the weak-star closure of E satisfying (3.2). In particular ||0||* ^ 1. 
Moreover, p(T) ^ lim^oo k~ l (j)(T (x) — x) ^ ||0||*p(T), and since p(T) = p(T) > 0, 
1 < ||0||*, which shows that ||0||* = 1. □ 

Actually, the assumption that p(T) > can be dispensed with in the finite dimension 

case: 

Corollary 21. Let (X, || • ||) be a finite dimensional Banach space, and let T : X — >■ X 
be non-expansive. Then for any x G X, there exists a linear form cf) of norm one, in the 
closure of the set of extreme points of the unit ball, such that 



<p(T k {x)) ^ 4>{x) + kp(T) 



(3.3) 



for all k G N. 



Proof. The proof is similar to the one of Corollary [20} We still take E as the set of extreme 
points of the dual unit ball, which belongs to the dual unit sphere, and we conclude since 
the unit sphere of a finite dimensional Banach space is closed. □ 

Remark 22. The case in which p(T) = 0, in finite dimension, has been studied in particular 
who proved that either T has a fixed point or there exists a linear form 
1 such that ip(T k (x)) — > +oo for all x G X. This result is established by 



by Lins in Lin09 
if of norm 

showing first that in the latter case, there is a horofunction h such that h(T k (x)) — > +oo 
as k tends to oo (modulo a change of sign convention for horofunctions). As noted in 
Remark 3.1 of |Lin09| , the construction of the linear form if from the horofunction h 
relies essentially on a sub-differentiability argument. An argument of the same nature is 



used here to derive Corollary 21 from Theorem 16 Note however that Corollary 21 and 



the existence of the horofunction in Lin09 are independent results (none of them can be 
recovered from the other one). 



3.3 Application to Shapley operators 

In this section Q is a compact topological space, and X = ^(Q) is the Banach space of 
continuous functions from Q to IR, endowed with the sup norm. We will be interested in 
functions T : X — > X that are order-preserving, meaning that 

x ^ y => T(x) ^ T(y) 
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where ^ denotes the canonical partial order on functions. We shall say that T commutes 
with the addition of a constant if 



T(A + x) = A + T(x), 



VA g 



where A + x denotes the function + x(oj). 

As in Example [9j we set t(x) := max ue j]i(u) and consider S(x,y) = t(y — x). Recall 
that S is a hemi-metric compatible with the sup norm and that (X, 5) is metrically star- 
shaped. 

It is easy to see that T is order-preserving and commutes with the addition of a 
constant if and only if 

t(T(x) - T{y)) <: t(x - y) 
for all x and y, meaning that T is non-expansive for S (see in particular GK95 , or LS05| ). 



For this special hemi-metric 5, we denote by p + (T) and p+(T) the quantities defined 



in (2.4) and (2.5). Theorem 16 and Corollary 18 imply that 



p + (T) = lim max 

fc— >+oo wGf! 



T k {x){uj) -x(u) 
k 



P+( T ) = in £ max(T(y)(u;) - y{u)) 



for all x in X. 



We next derive from Theorem 19 the following generalisation of a result established 
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GG04 when Q is finite. 



Corollary 23. Let Q be a compact topological space, and let X = ^(Q) be the Banach 
space of continuous functions from Q to IR, endowed with the sup norm. Suppose that 
T : X — > X is order-preserving and commutes with the addition of a constant. Then for 
every x G X , there exists uj + 6 Q such that, 



T\x){uj + )^x{uj + ) + kp + {T) 



(3.4) 



for all k E N. 



Proof. We apply Theorem 19 to the hemi-norm p(z) = t(z), taking for E the set of 
evaluation functions e u , 00 G Q. Since Q is compact, the set E is weak-star closed, so that 
the linear form <fi in Theorem 19 can be written as = e u+ for some u + G E. □ 



Order-preserving maps T commuting with the addition of a constant arise as Shapley 
operators in game theory. Suppose now that Q, A and B are three compact metric spaces, 
that g is a continuous function from A x B x Q to M, and that q : A x B x Q — > A(fl) is 
continuous, where A(f2), the set of probabilities on Q, is endowed with the weak topology. 
Then the Shapley operator of the two-player, zero-sum stochastic game with state space 
Q, action sets A and B, payoff g, and transition probability q, is defined as 



T(x)(lo) = sup inf < / 

M eA(A) "6A(B) {J AxB 

= inf sup < / 

veA(B) M gA(4) {JAxB 



g(a,b,oj) + / q(dw'\a,b,uj)x(uj') 



g{a,b,uS) + / q(dw'\a, b, u)x(u f ) 
Jn 



p(da)v(db) 
n(da)v{db) \ . 
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Th e map T is order -preserving and commutes with the addition of a constant. It is shown 



in 



MSZ94 



Now03 



that it preserves ^(O). Then, T represents the average reward 

per time unit if both players play optimally in the fc-stage game and if the starting state 
is oj. Thus, Corollary [23] shows that there is an initial state u + which is the best for the 



maximising player in the long term. In particular 



converges as k goes to infinity. 



Remark 24. To every result concerning the hemi-norm t corresponds a dual result con- 
cerning the hemi-norm b(x) := —t(—x) = min we n x(oS). In particular, 



lim min 

k— >+oo u}£Q 



T k (x)(u)-x(u) 
k 



p_{T) := supmin(T(y)(o;) - y(u)) 



Moreover, there exists w_ G Q satisfying 

T fe (x)(w_) x{uJ) + kp_(T) 



for all fceN. This is readily deduced by applying Corollary 23 to the map x -» —T(—x). 



3.4 Comparison with a theorem of Karlsson 



Theorem 16 should be compared with the following result. (Recall that a metric space is 
proper if every closed ball is compact.) 



Theorem 25 (Karlsson |Kar01 , Th. 3.3]). If T is a non-expansive self-map of a proper 
metric space (X, 5), then, for all x G X, there exists a horofunction h (depending on x), 
such that 

h(T k (x)) ^ h{x) + kp(T) , 

holds for all k ^ 1 . 



Note first two differences in the statements of Theorem [16] and of Karlsson's theorem: 
we require a metric convexity assumption whereas Karlsson does not; Karlsson makes 
the general assumption that the metric space is proper whereas we do not. Actually, the 
compactness issue appears to be a mere technicality: it can be checked that the proof 
of Theorem [25] remains valid even if X is not proper if one defines horofunctions with 



respect to the topology of pointwise convergence as we do here (see Remark 14). 



A key difference however lies in the conclusions of both results. In Theorem 16, horo- 
functions are independent of the choice of x, whereas in Karlsson's result, the horofunction 



does depend on x. Note also that inequalities involving piT) as in Theorem 16 and Corol- 
lary 17 are stronger than their counterparts involving p(T) since p(T) ^ p{T). This is 



illustrated in Example [26] below, in which we consider a non-star-shaped space for which 
Theorem [25] applies but not Theorem [16} This example shows that the assumptions of 



Theorem 16 cannot be weakened 
shaped. 



it is primordial that the space X is metrically star- 



Hence, Theorem 16 and Karlsson's result are incomparable: the latter holds 
under more general circumstances, but the former yields a stronger conclusion. 

Example 26. Let T := IR(modl) be the torus with its canonical metric dj; and let X = 
Txi, with the metric 



d((x, t), (x', t')) := dj(x, x) + \t — t'\ 
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Let a G (0, |) \ Q, and consider the mapping T : X — > X, 

T(x, t) — (x + a, t + 1) 
which is non-expansive for the metric d. It is straightforward to check that p(T) = 1 



p(T) = p(T) by Corollary 18). 



but p(T) = 1 + a, so inequation (3.1) cannot be satisfied (otherwise one would have 



In fact we prove that it is not even true that there exists a function h G ^# such that 
for all x, 

h(T k (x, t)) ^ h(x, t) + kp{T), \/k G N. (3.5) 

Fix a basepoint (x,t) in X. Since the metric on X = T x R is the sum of the two 
metrics (T, dy) and (R, | • |), we can write any Martin function /i of X as the sum 

h(x,t) = hi(x) + /i 2 (^) 

where hi G is a Martin function of (T, dx) and h 2 is a Martin function of (R, | ■ |), 



relatively to the basepoints x and i, respectively (see Prop. 9.1 in AGW09) ). The only 
Martin functions on T are hi(x) = —dj(x, y) + df(x, y) with y G T. The Martin functions 
of (R, | ■ |) are either of the form h 2 (t) = —\t — s\ + \i — s\, with s G R, or the affine 
functions h2(t) = t — i and h 2 = —(t — i). We search a Martin function h such that, for 
any (x, t) G X, 

h(T k (x, t)) ^ h(x, t) + k, Vk G N (3.6) 



i.e., 



h 1 (x + ka)+h 2 {t + k)^h 1 {x) + h 2 {t) + k, Mk G N. (3.7) 

Any Martin function h\ is bounded. If h 2 is of the form h 2 (t) = —\t — s\ + \t — s\ or 
^2^) — — (t — t), then, for a fixed (x,t), the left-hand side of (3.7) remains bounded 
from above as k — > oo, whereas the right-hand side of (3.7) tends to infinity. The only 
possibility is thus h 2 (t) = t — i. In that case, ( 3.7[ ) becomes 



hi(x + ka) — hi(x) ^ 0. 

Since the sequence (x + ka)ken is dense in T, and since hi is continuous, we have that 
^i(^) ^ inf^gT hi(y). This is true for any x G T, so hi is constant. But there are no 
constant Martin function of (T, dj). 

Karlsson's result shows however that for any choice of (x,t) G X, there exists a 



horofunction h (depending on (x, t)) such that (3.6) is satisfied 



4 A Denjoy- Wolff type theorem for order-preserving 
homogeneous self-maps of symmetric cones 

4.1 The reverse Funk metric 

In this section, we recall or establish preliminary results concerning cones. We consider a 
(closed, convex and pointed) cone C in IR n with non-empty interior X := int C . We equip 
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C with the partial order defined by x ^ y if y — x G C. For all x, y G C \ {0}, we set 

M(y/x) := inf{A > | Ax ^ y}, m(y/x) := sup{A > | Ax ^ y} . 

Since C is closed, as soon as the two latter sets are non-empty, their respective infimum 
and supremum are attained; in particular, M(y/x) and m(y/x) are finite and positive. 
The (reverse) Funk hemi-metric on the interior of C is defined by 

RFunk(x, y) := logM(y/x) . 

More generally, we shall use the notation RFunk(x, y) as soon as x, y G C \ {0} are such 



that fix ^ y for some \i > 0. The terminology "Funk metric" is borrowed from PT09 



Wal08 ; it refers to |Fun29 . The map 8 := RFunk is easily checked to be a hemi-metric on 



int C. Indeed: 5(x, y) is finite since x is in the interior of C; the triangular inequality and 
the fact that 5(x, x) = for any x G X are easy to verify since C is pointed; the fact that 
S(y,x) = S(y,x) = ==>■ x = y follows from the fact that C is closed and pointed. The 
term "reverse" is by opposition to the hemi-metric (x, y) t— > 5(y, x) which has a different 



horoboundary Wal08|. 



The quantities M(y/x) and m(y/x) can be expressed in terms of the extreme rays of 
the dual cone C* . The latter is the set of continuous linear forms taking non-negative 
values on C . For each cone K, Extr K denotes the set of representatives of the extreme 
rays of K (i.e., the non-zero vectors belonging to these extreme rays). 

Lemma 27. Let x,y,z G C \ {0} be such that \ix ^ y and z ^ vx for some //, v > 0. 
Then, 

M(y/ X ) = sup sup (4.1) 

m{z/x) = inf ^4 = inf 44 ■ (4-2) 

<^eC*, <p(x)^0 <f(X) i^eExtrC*, <p(x)^0 (p(x) 

Moreover, when x G int C, the latter suprema are attained, and the condition that (p(x) ^ 
can be replaced by tp G C* \ {0}. When z G intC, the latter infima are attained. 

This result is somehow standard. We include the proof for the convenience of the 
reader. 

Proof. Since C = (C*)*, the inequality Ax ^ y, i.e., Ax — y G C, is equivalent to <^(Ax — 



y) ^ 0, for all <p G C*, from which the first equality in (4.1) follows. 

Let now recall that if u is in the interior of C, and if (p G C* \ {0}, then, <p(u) cannot 
vanish. Indeed, we can find a ball B centred at 0, such that u + z G C for all z G B. If 
<p(u) = 0, then, ^ ip(u ± z) = ±<p(z) holds for all z G -B, and so = 0. 

Choose now a vector w G intC, and define S u := {(p G C* | = 1}. Observe 

that S u is compact and convex, and that each extreme ray of C can be identified to the 
unique extreme point of E u which generates this ray. Let A denote the value of the last 



supremum in (4.1), so that the inequality \<p(x) is (p{y) holds for all extreme points ip of 



S u . Since every element (p of S n is a barycentre of extreme points of S u , we deduce that 
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the same inequality holds for all points ip of E U) which shows that A coincides with the 



first supremum in (4.1). 



The arguments to establish (4.2) are dual 



Finally, when x G intC, we consider the map J : Yi x — > R, <p i— > <p(y)/(p(x) = (f(y)- 
Since J is linear, it attains its maximum at an extreme point of the compact convex set 
£ x , which implies that both suprema in (4.1) are attained. Similarly, when z G intC, we 
note that every ip G S z maximising attains the first infimum in (4.2), and so, the 
second infimum is attained at an extreme point of T, z . □ 



Walsh showed in Wal08 that the horoboundary of int C in the (reverse) Funk hemi- 
metric coincides with the space of rays contained in the Euclidean boundary of C. 



Proposition 28 ( |Wal08[ Prop. 2.5]). Let C be a closed convex pointed cone of non- 
empty interior in a finite dimensional Banach space, with basepoint x. Then, any Martin 
function h for the (reverse) RFunk hemi-metric corresponds to a vector u G C \ {0} ; 



h(x) 



RFunk(x, u) + RFunk(a;, u) ,Wx £ int C, 



(4.3) 



and h is a horofunction if and only if u G dC \ {0}. 



We shall also be interested in the following variant of the Funk hemi-metric, considered 
by Papadopoulos and Troyanov |PT09) : 

RFunk" 1 " (x, y) := max(RFunk(x, y), 0) . 

The map RFunk" 1 " is easily seen to be a hemi-metric. The following simple application of 
the arguments of Wal08 characterises the boundary of the hemi-metric RFunk + . 



Proposition 29. Let C be a closed convex pointed cone of non-empty interior in a finite 
dimensional Banach space. Then, the Martin space for the RFunk" 1 " hemi-metric consists 



of the functions (4.3), with u G C \ {0}, together with the functions 

h(x) = — RFunk + (x, u) + RFunk + (x, u) , Vx G int C, 



(4.4) 



with u G C , and an element is a horofunction if and only if it is of the form (4.3), with 



u G C \ {0} ; or of the form (4.4) with u G dC 



Proof. Any function h of the Martin space is the pointwise limit of a sequence of functions 
h m '. x i — y — RFunk + (x, z m ) + RFunk + (x, z m ), where (z m ) m ^i is a sequence of elements 
of intC. If \\z m \\ — > oo as m tends to infinity, then, for all x G intC, RFunk ( x, z m ) also 
tends to infinity, and so RFunk" 1 " (x, z m ) = RFunk ( x, z m ) for m large enough. Then, the 
result of |Wal08| shows that h = lim m h m is of the form (4.3), for some u G C\ {0}. It 



remains to examine the case in which z m contains a bounded subsequence. Then, possibly 
after replacing z m by a subsequence, we may assume that z m converges to a point u G C. 
Using the fact established in Wal08| that for x G int C, RFunk(x, y) is continuous in the 
second argument y G C, we deduce that h m converges to the function (4.4). Finally, the 
characterisation of the horofunctions is straightforward. □ 
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The metric dx associated to the Funk hemi- metric, 



cIt{x, y) = max(RFunk(x, y), RFunk(y, x)) 



is called the Thompson metric. We recall that for this metric, the map 

(3 s - a s (3 s - a s 



(3a s - a/3 s + (3 s - a 



-x" + 



(3a s - a(3 s + (3 s - a 



(4.5) 



with (3 = exp(RFunk(x°, y)) and a = exp(— RFunk(?/, a; )) is a (straight line) geodesic 
joining x° to y, see |Nus88| . However, as shown in [N W04] , this choice of geodesies does not 
satisfy Busemann's non-positive curvature condition (even when C = IR+). Actually, it is 



shown in NW04 that only a weaker inequality is valid. The latter inequality is optimal 



(it is straightforward to find explicit counter examples to the convexity inequality). 

However, when C is the cone S„ of n x n positive semi-definite matrices, there is 
a different choice of geodesies between any pair of matrices Z, Y in the interior of C, 
involving geometric means: 



1y{s) 



z~ l 2 ( z-Wz- 1 ? 



z\ 



(4.6) 



See in particular Nus88 , NW04 . The latter are indeed known to be geodesies in Thomp- 



son's and Hilbert's metric. Note that when Z is the identity matrix, this geodesic seen 
with logarithmic glasses is a straight line: log7y(s) = slogF. The general form of the 
geodesic is obtained by considering the linear transformation X i— > Z~ 1 ^ 2 XZ~ 1 ^ 2 , which is 
an automorphism of the cone of positive definite matrices, isometric both in Thompson's 
and Hilbert's metric, and sending Z to the identity. This construction actually makes 
sense more generally when C is a symmetric cone [FK94] , meaning; that C is self-dual 
{C = C*), and that the group of automorphisms leaving invariant the interior of C acts 
transitively on this interior. Symmetric cones include in particular the Lorentz cones, the 
cones of real or complex positive semidefinite matrices, and Cartesian products of such 



cones. It is shown in GW03 that the metric convexity property of Definition 5 holds for 



any symmetric cone, for the latter choice of geodesies, for any choice of the centre, both 
for the Thompson and Hilbert metric. A similar observation was made by Lawson and 



Lim LL07b . Actually, the arguments in GW03| apply as well to the Funk and RFunk 



hemi-metrics. Hence, the present results apply to those self-maps of the interior of a 
symmetric cone that are non-expansive in any of the previously mentioned hemi-metrics. 

Remark 30. When C = S£ is the cone of positive semi-definite matrices, when x G intC, 
and y G C\ {0}, M(y/x) coincides with the maximal eigenvalue of the symmetric matrix 
v := x~ l l 2 yx~ l l 2 . Moreover, denoting by w the (rank one) orthogonal projector on any 
eigenline corresponding to this eigenvalue, one can check that the linear form ip(u) := 
tr(x~ 1 / 2 wx~ 1 / 2 u) attains the supremum in (4.1). When x G C \ {0}, with y ^ [ix for 



some /i > 0, a maximising linear form ip can be constructed in the same way, replacing 
by the square root of the Moore-Penrose inverse of x. A dual argument allows one to 



x 



-1/2 



construct a linear form attaining the infimum in (4.2). Note also that these constructions 



can be easily adapted to the case of a symmetric cone C. 
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Remark 31. In the case of the cone of symmetric positive definite matrices, the Thompson, 
Hilbert, Funk, and RFunk + hemi-metrics all have the following form 

8„(A, B) := v{\og Spec(A- 1 5)) = p{\og A a , . . . , log A n ) A, B e int S+ , 

where v is a symmetric function IR n — y R (by symmetric, we mean that v is invariant by 
permutation of its variables), Spec(M) denotes the sequence of eigenvalues of a matrix 
M, counted with multiplicities, so that Ai, • • • , A n denote the eigenvalues of the matrix 
A~ 1 B (which are real and positive since A and B are positive definite). The previously 
listed hemi-metrics are obtained by taking for v the following maps, W 1 — y R, 



a i — y max uJ 



/i i — y max \fii — jjLj\ 



ji I—)- max /ij, 



H i — y max(max /ij, 

l<i<n 



0) 



respectively. We may consider more generally any symmetric hemi-norm v. Then, the 



corresponding map 5 U can be checked to be a hemi-metric, which is of Finsler type Nus94 



Indeed, measuring the length of a vector Y in the tangent space of S£ at point X by the 
hemi-norm u(Spec(X~ 1 Y)) , we get 



S V (A,B) 



inf 

P 



z/(Spec(/3(s)- 1 /3(s))c/s 



the infimum being taken over the set of curves /3 from [0, 1] to the interior of the cone, such 
that /3(0) = A and /3(1) = B. We recover as special cases of this construction the invariant 
Riemannian metric |Mos55| {v is the Euclidean norm), the log p-Schatten metric |FF04| 



[y is the £ p norm), as well as metrics arising from symmetric Gauge functions Bha03 (the 



latter are symmetric hemi-norms that are invariant by a change of sign of each variable, 
and take positive values except at the origin). For any choice of symmetric hemi-metric 
v, the map (4.6) yields a family of geodesies for the hemi-metric 5„ joining a centre Z 



to any matrix Y, for which int S£ is a metrically star shaped space. This can be proved 
along the lines of |Bha03| , where the same property is established in the special case in 
which v is a symmetric Gauge function. The reader may also consult ACS00 GW03I, in 



which the significance of several matrix inequalities in terms of non-positive curvature is 
brought to light. 



A number of non-linear maps, including discrete Riccati operators (Example 40 below), 
are known to be non-expansive in various metrics 5 Ul see 



Bou95||LW94l|LL08HLL07al 



4.2 A Denjoy- Wolff type theorem 

We shall now consider specifically the case of order-preserving and positively homogeneous 
maps. A self-map T of C is order-preserving if x ^ y T(x) ^ T(y). It is positively 

homogeneous (of degree one) if T commutes with the product by a positive constant. 
These notions make sense more generally if T is only defined on a subset of C invariant 
by any such product. It is readily checked that a self-map T of the interior of C is order- 
preserving and positively homogeneous if and only if it is non-expansive in the Funk 
metric. 
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Corollary 32. Let T be an order-preserving and positively homogeneous self-map of the 
interior of a symmetric cone C , and let 



p(T):= inf RFunk(y, T(yj) . 

y£mt C 

Then, for all x G int C , 

- P (T) = p(T) := lim 

fc— >oo k 

and there exists an extreme ray of C , such that for any element w in this ray, 

log(w,T k (x)) ^ \og(w,x) + kp(T). 
Here, (•, •) denotes the scalar product of the Euclidean space in which C is self-dual. 



Proof. The map T is non-expansive in the Funk metric. By Corollary 18, p(T) = p(T). 



Let h be the Martin function appearing in Theorem 16 Then, there exists a sequence of 
elements of int C, (z m ) m ^i such that i(z m ) converges to h. Let us now fix x G int C . Since 
C = C*, any element (p of an extreme ray of C* can be written as tp(y) = (w,y) where 
w is an element of an extreme ray of C. Then, it follows from Lemma [27] that for all m, 
there exists an element w m G Extr C such that 

RFunk(x, z m ) = log -. 

{w m ,x) 

A ray of C is known to be extreme if and only if it contains a primitive idempotent of 



the associated Jordan algebra FK94, Prop. IV.3.2]. Moreover, the set of these primitive 
idempotents is compact (Exercise 5 of Chap. 4 in |FK94| ). Since the previous property 
is invariant by a scaling of w m , we may assume that each w m is a primitive idempotent, 
and so, possibly after extracting a subsequence, we may assume that w m converges to a 
primitive idempotent w, which therefore belongs to an extreme ray of C. For all k ^ 0, 
we get 

log(w,T fc (x)) - log(w,x) = limlog(w m ,T fc (x)) - log(w m ,x) 

m 

^ lim - RFunk(T fc (x), z m ) + RFunk(x, z m ) 

m 

= h{T k {x)) - h{x) ^ kp(T) 
which concludes the proof. □ 



We next point out a variant of Corollary [32] concerning order-preserving and sub- 
homogeneous self-maps of the interior of a symmetric cone. Recall that sub-homogeneous 
means that for all x G intC and A ^ 1, T(Xx) ^ XT(x). Order-preserving and sub- 
homogeneous maps are easily seen to be characterised by the following property 

RPunk + (T(x),T(y)) < RFunk + (x,?/) . 
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Corollary 33. LetT be an order-preserving and sub-homogeneous self-map of the interior 
of a symmetric cone C , and let 



p(T) := inf RFunk+(y,T(y)) 

ySint C 



Then, for all x G int C , 



lim 

k— >oo 



RFunk + (a;,T fc (a;)) 
k 



and, if p(T) > 0, there exists an extreme ray of C, such that for any element w in this 
ray, 

\og(w,T k (x)) ^ \og(w,x) + kp(T). 
Proof. Since the hemi-metric RFunk + is bounded from below, the observation made at 



the end of the proof of Theorem 16 shows that if z m is any sequence of points of int C such 
that i(z m ) converges to the Martin function h in this theorem, then, the sequence z m must 
be such that ||-z m || — > oo. Then, it follows from the characterisation of the boundary of 



the RFunk metric (proof of Proposition 29) that h, which cannot be of the form (4.4), is 
necessarily of the form (4.3). Therefore, we conclude as in the proof of Corollary 32 □ 



Remark 34. Corollary [32] should be compared with a related result of Gunawardena 
and Walsh |GW03[ Th. 1], in which the map T is only required to be non-expansive 
in Thompson's metric. Then, it is shown that all the accumulation points of the sequence 
{\ogT k (x))/k belong to the same face of a ball which respect to the norm in the tangent 
space of C at point x. 



Remark 35. Lins established a Denjoy- Wolff theorem of a different nature in Lin07 



concerning self-maps of a polyhedral domain that are non-expansive in Hilbert's metric. 



4.3 The Collatz-Wielandt theorem 



We now show that the Collatz-Wielandt formula (1.2) is a special case of the maximin 



formula (1.1). 



If T is an order preserving and positively sub-homogeneous self-map of the interior of 
a (closed, convex and pointed) cone C in a finite dimensional vector space, we define the 
radial extension T of T to the closed cone C by 

fix) := lim T(x + ez), Vx G C , 

e->0+ 



where z is any point in the interior of C. One verifies that this definition is independent 
of the choice of z\ equivalently, we may define T(x) as mi{T(y) \ y — x G intC}. The 
order preserving and positively (sub-) homogeneous character of T carries over to T . It is 



shown in BNS03 1 that T does not necessarily have a continuous extension to C, but that 



it always does if C is polyhedral. 
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Lemma 36. Let T be an order-preserving and positively sub-homogeneous self-map of the 
interior of a (convex closed and pointed) cone C in a finite dimensional vector space, and 
consider the Martin function h in (4.3) associated to a vector u G C\ {0}. Then, for all 
AgI, 

(Vx G intC, h{T{x)) - h(x) ^ A) <=> (V 7 > 0, f ( 7 u) ^ 7exp(A)w) . (4.7) 

Proof. Since, h(T(x)) — h(x) = — RFunk(T(x), u) + RFunk(x, w), the left-hand side con- 
dition in (4.7) holds if and only if 

u < exp ( - A + RFunk(x,w))T(x), Vx G int C . (4.8) 

For e > and 7 > 0, define x e := 7W + ez, where z is a given element in int C. Then, 
RFunk(x e ,u) ^ RFunk(7w, m) = — log 7. Taking x = x e in (4.8), and letting e tend to 
zero, we get 

u ^ exp(— A)7 _1 T(7m) , 



and so, the condition at the right-hand side of (4.7) holds. Conversely, fix x G intC, and 
choose 7 such that 7RFunk(x,w) ^ 1. We have 



7U ^ exp(— X)T(ju) (using the right-hand side of (4.7)) 

^ exp(— A)T(7exp(RFunk(x, u))x) (since T is order preserving) 
^ exp(— A)7 exp(RFunk(x, u))T(x) (since T is sub-homogeneous) 
= 7exp(— A + RFunk(x, u))T(x) . 

Cancelling 7, we arrive at ( |4.8 ). □ 

Observe that when T is positively sub-homogeneous, the following recession map 

f r :C^C, f r (x) = lim 7 _1 T( 7 x) = inf 7 _1 T( 7 x) (4.9) 

7— >oo 7>0 

is well defined. Of course, T = T r if T is positively homogeneous. We now arrive at the 
following corollary of Theorem [TJ which extends the Collatz-Wielandt formula (1.2). 

Corollary 37. Let T be an order-preserving and positively sub-homogeneous self-map of 
the interior of a symmetric cone in a finite dimensional vector space, and assume that 

;= Um logRFuukferV)) > Q 

Then, 

■ t , (w,T(y)) . (w,f r (u)) 

ml max log — r — = p(l ) = max mm log — r — . (4.11) 

yGintC MJSExtrC (w, y) «eC\{0} weExtrC* \W,U) 

Moreover, when T is positively homogeneous, the same conclusion remains valid even 
when the condition p(T) > does not hold. 
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Proof. We apply Theorem 16 to the space X := intC equipped with the hemi-metric 
5 := RFunk + . Let us denote by p+{T) the escape rate of T in this hemi-metric, so that 



P+ (T)= inf RFunk + (y, T(y)) . 



(4.12) 



Since RFunk + = max(RFunk, 0), the hemi-metrics RFunk and RFunk + coincide as soon 
as RFunk + or RFunk is positive. Hence, p(T) > implies that p+(T) = p{T) and that the 



term RFunk (y,T(y)) in (4.12) can be replaced by RFunk(y, T(y)). Together with (4.1) 
this gives the first equality in (4.11). 



We observed in the proof of Corollary 33 that when the escape rate in the RFunk 



hemi-metric is positive, the corresponding horofunction h in Theorem 16 is necessarily of 



the form (4.3). Since the limit coincides with the infimum in the definition of the recession 



function (see (4.9)), the right-hand side condition in (4.7) is equivalent to 



T r (u) ^ exp(A)ii 



Then, using Lemma 36 , we rewrite the maximin characterisation in Theorem 16 



as 



p(T) = max inf h(T(x)) — h(x) 

h x&ntC 

= max{A | 3u G C \ {0}, T r (u) ^ exp(A)w) 

(w,f r (u)) 



max log inf 

«eC\{0} weExtrC (w,u)^0 \W,U) 



(by (4.2)). 



Moreover, it follows from Remark 30 that the latter infimum is attained. 



Finally, when T is positively homogeneous, we equip the space int C with the hemi- 



metric RFunk instead of RFunk . Then, the first equality in (4.11) follows readily from 



Theorem 16 and (4.1) (the assumption that piT) > is not needed any more), and the 
proof of the second inequality is unchanged. □ 



Remark 38. The characterisation (4.11) can be rewritten equivalently as 



p{T) = log inf {p > | 3y E mtC, T{y) ^ py} 

= logmaxjp ^0\3ueC \ {0}, T r {u) ^ pu} . 



(4.13) 
(4.14) 



Remark 39. When C is the standard positive cone, and T is a continuous order preserving 
and positively homogeneous self-map of the closed cone C, Nussbaum |Nus86| showed that 

p(T) = logmaxjyU ^ | 3u G C \ {0}, T{u) = pu} , 



which is more accurate than (4.14). Moreover, some of the results of |Nus86| have been 



extended in |MPN02| to more general (normal, possibly infinite dimensional) cones. How- 
ever, Corollary 37 is applicable when T is only defined on the interior of the cone and 
does not extend continuously to its boundary (this situation does occur in practice, see 
the discussion in |BNS03] of the matrix harmonic mean). 
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Example 40. We finally illustrate the previous Denjoy- Wolff type results with an example 
originating from quadratic optimal control. Let A, B be n x n positive semidefinite matri- 
ces, and let M be a n x n invertible matrix. Consider the following discrete time Riccati 
operator, which is a self-map of the interior of the cone S£ of n x n positive semidefinite 
symmetric matrices, 

T(X) = A + M(B + X- l y l M* , 
where M* denotes the transpose of M. We refer the reader to Bou95 , LW94 , LL08I for 



more background on these maps. In particular, T is order-preserving, and, as shown by 



Liverani and Wojtkowski LW94|, T is non-expansive in the Thompson metric. We noted 



above that the latter property, together with the preservation of order, implies that T is 



positively sub-homogeneous, and so, Corollaries 33 and 37 apply to the map T. 

If the matrices A and B are positive definite, then T is known to be a strict contraction 



in Thompson metric LW94 1 . It follows that T has a unique fixed point in int S£ to which 



every orbit converges. In particular, the escape rate p(T) is zero. 

We next examine a more interesting case in which B is of rank one, leading to a 
positive escape rate. Assume that n = 2, write B = vv* for some (column) vector v, and 
assume that M is of the form M = al, with a > 1, where / is the identity matrix. 

We claim that p(T) = 2 log a. To see this, let u denote a non-zero vector orthogonal to 
v, and consider U = uu* . Then, it is easily checked that the radial extension of T satisfies 
f(jU) = A + a 2 ^U for all 7 > 0, and so, f r (U) = a 2 U. Using (|4.14[), we deduce that 



p(T) > 2 log a. To show that the equality holds, we observe that T(X) ^ A + MXM*. 
Hence, for all s > 0, T(sl) < A + a 2 sl < (Xi(A)s 



-1 



a 



)sl, where Xi(A) denotes 
the maximal eigenvalue of A, which by (4.13), implies that p(T) ^ log(Ai(y4)s _1 + a 2 ). 



Letting s tend to infinity, we arrive at p(T) ^ 2 log a, and so p(T) = 2 log a. 



A possible choice of the horofunction h appearing in Theorem 16 



is 



h(X) = - RFunk(A, U) + RFunk(/, U) , 
where the basepoint is the identity matrix I. Then, the linear form x 1— > (w, x) constructed 



in the proof of Corollary 33 can be checked to be X 1— )■ (U,X) (or any scalar multiple of 
this form). 

For n = 2, the cone of positive semidefinite matrices has the shape of a Lorentz 
(ice-cream) cone. Then, the superlevel set h(X) ^ A is nothing but the interior of the 
translated Lorentz cone exp(A — RFunk(J, U))U + S 1 ^. This is illustrated in Figure [l] 
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